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We present numerically exact predictions of the periodic and single-impurity Anderson models to 
address photoemission experiments on heavy Fermion systems. Unlike the single impurity model 
the lattice model is able to account for the enhanced intensity, dispersion, and apparent weak 
temperature dependence of the Kondo resonant peak seen in recent controversial photoemission 
experiments. We present a consistent interpretation of these results as a crossover from the impurity 
regime to an effective Hubbard model regime described by Nozieres. 



Introduction Metallic compounds containing rare 
earth elements with partially filled / shells, such as 
CeBei3 or UPt3, belong to the general category of heavy 
fermions jjj. They are characterized by a large Pauli 
susceptibility and specific heat as compared to ordi- 
nary metals, which indicate a huge effective electronic 
mass, and also by anomalous transport properties such 
as non-monotonic temperature dependence of the resis- 
tivity. These anomalies are usually attributed to the for- 
mation of a resonant state at the Fermi energy due to 
the admixture of the electronically active and highly lo- 
cal / orbitals with the metallic band of the host. Heavy 
fermions are usually modeled by the single impurity An- 
derson model or the periodic Anderson model depending 
on the concentration of the correlated / orbitals. 

Photoemission experiments provide a direct probe of 
the characteristic resonant states in these materials. 
However, there are two apparently contradictory sets of 
photoemission results. The large body of spectroscopic 
data accumulated by the Bell Lab and Neuchatel groups 
for scraped polycrystalline samples consistently indicate 
qualitative and sometimes quantitative agreement with 
the predictions of the single impurity models Q . On the 
other hand some recent photoemission results obtained 
by the Los Alamos group for single-crystalline samples 
cleaved in situ have revealed several apparent inconsis- 
tencies with the predictions of single impurity models, 
including enhanced intensity, a much weaker tempera- 
ture dependence and a dispersion of the Kondo peak || . 
This has caused a persistent dispute regarding the ability 
of single impurity models to describe the observed reso- 
nant state in Kondo materials Q . We refer the reader to 
Malterre et al. ^ for a recent review. 

The single impurity Anderson model (SIAM) is typi- 
cally used as the paradigm for comparison with exper- 
iments due to its universality, the abundance of good 
approximations, and the existence of exact solutions. Al- 
though at high temperatures the SIAM captures the same 
physics as the lattice model, it cannot account for the 
electronic coherence at low temperatures. The periodic 



Anderson model (PAM) is believed to correctly describe 
the strong correlation of / electrons as well as their co- 
herence at low temperatures and the interaction between 
the screened moments. 

We have recently analyzed some static properties of the 
PAM §. We found that when the / band filling n f » 1 
and the d band filling nd ^ 0.8 (so that the system is 
metallic), the Kondo scale for the PAM, To, is strongly 
suppressed compared to Tq IAM , the Kondo scale for a 
SIAM with the same model parameters. However the 
high temperature (T > Tq IAM ) properties of the two 
models are similar, so that Tq IAM is also the relevant 
scale for the onset of screening in the PAM whereas T is 
the scale for the onset of coherence. We found no univer- 
sal relation between Tjf IAM and T (their ratio depends 
on the d band filling) . We also demonstrated that in the 
screening regime (T < Tq IAM ) the rate of change of the 
local magnetic moment with temperature is smaller in 
the PAM than in the SIAM. 

In this paper we show that the PAM predicts a much 
weaker temperature dependence for the Kondo peak than 
the SIAM, consistent with the experiments on single- 
crystalline Kondo lattice materials. We also show that 
the Kondo peak is dispersive in the PAM, giving rise to 
heavy quasiparticle bands near the Fermi energy. We 
use Nozieres' idea of the effective Hubbard model for the 
screening clouds to interpret this "band formation" 
and the slow evolution of the Kondo peak. This also 
gives insight into the emergence of two relevant energy 
scales (Tq IAM and T ) in the PAM. 

Method The PAM Hamiltonian on a D-dimensional 
hypercubic lattice is, 
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In Eq. (1), d(f)\J destroys (creates) a d(f) electron with 
spin a on site i. The hopping is restricted to the nearest 
neighbors and scaled as t = t* /2\fT>. U is the screened 
on-site Coulomb repulsion for the localized / states and 

V is the hybridization between d and / states. This 
model retains the features of the impurity problem, in- 
cluding moment formation and screening, but is further 
complicated by the lattice effects. 

Metzner and Vollhardt || observed that the irreducible 
self-energy and vertex-functions become purely local as 
the coordination number of the lattice increases. As a 
consequence, the solution of an interacting lattice model 
in D = oo may be mapped onto the solution of a local 
correlated impurity coupled to a self-consistently deter- 
mined host ||. We employ the quantum Monte Carlo 
(QMC) algorithm of Hirsch and Fye to solve the re- 
maining impurity problem and calculate the imaginary 
time local Green's functions. We then use the maximum 
entropy method |fl] to find the / and d density of states 
and the self-energy Q . 

Results We simulated the PAM for a wide variety of 
fillings and parameters in units of t* (considered to be a 
few electron-volts, the typical band-width of conduction 
electrons in metals). Here we present the results for U — 
1.5, V — 0.6, 71/ 1 with two different d band fillings 
iid = 0.4 and rid = 0.6 for which the Kondo scales are 
T = 0.014 and T = 0.054 respectively. We use the 
symmetric limit of the SI AM (n/ = rid = 1) for the 
comparison since the results for the SIAM are universal 
and hence independent of the filling. 

Fig. 1 shows the / density of states for the PAM and 
SIAM near the Fermi energy Ep = 0. Both models show 
temperature-dependent Kondo peaks of similar width, 
but the rate of change of the peak intensity with tem- 
perature is much smaller and the Kondo peak persists 
up to much higher temperatures in the PAM than for 
the SIAM. This is consistent with the protracted behav- 
ior of screened moments in the PAM ||, i.e., the rate 
of change of the screened moments with temperature is 
much smaller in the PAM compared to the SIAM when 
rid ^ 0.8 and rif w 1. Note the difference in scales of 
the two parts of this figure, both the intensity and the 
spectral weight of the Kondo peak are larger in the PAM 
than in the SIAM although the hybridization parameter 

V is larger for the PAM in this case. This shows (and we 
generally find) that the height of the Kondo peak in the 
PAM does not scale like 1/V 2 as it does in the impurity 
models. For the PAM, starting at high temperatures the 
Kondo resonance peak is located well above the Fermi 
energy but as we lower the temperature it shifts toward 
Ep. Hence there is a spectral weight transfer into the 
region around the Fermi energy as the temperature is 
lowered but we can see that the spectral weight in this 
region increases slower in the PAM than in the SIAM. 

Fig. 2 shows the momentum dependence of the / and d 



spectral functions for the PAM along the diagonal direc- 
tion of the Brillouin zone (the main conclusions do not 
depend on the chosen direction) . In the large- D hypercu- 
bic lattice, the zone center (corner) corresponds to a very 
large negative (positive) e^. We see that near the zone 
center (lowest part of the figure) there are two apparent 
maxima, the lower one has mostly d character the upper 
one mostly / character. The latter has a narrow Kondo 
like feature but would not be seen in photoemission ex- 
periment since it is located above the Fermi energy. The 
quasiparticle peak (Kondo peak) near the Fermi surface 
starts to develop only as k moves well away from the 
zone center. Note that there is a gradual transformation 
of this peak from a mixed /d-character at e& = to an 
almost entirely / character at and above = 1.23. In 
addition to these narrow peaks the data shows a small 
and broad non-dispersive peak near ui = —0.8 which has 
mostly / character. This is a remnant of the lower un- 
hybridized / level, we were not able to resolve the peak 
corresponding to the upper / level for these model pa- 
rameters (we resolve this peak if we use a larger value for 
U/V 2 ). 

We clarify this situation in Fig. 3. The symbols in this 
figure show the position of / and d spectral functions 
maxima versus and the solid line shows the quasipar- 
ticle energy, calculated from the real part of the pole of 
the Green's functions. In the narrow region above the 
Fermi surface we see an almost dispersionless band. The 
imaginary part of these poles are very large, so they corre- 
spond to very broad peaks in the spectral functions which 
we were not able to resolve. The bands above and below 
this region correspond to well defined peaks in the spec- 
tral functions. The general features of this band struc- 
ture for the PAM persists up to very high temperatures 
(T/T > 10). 

The imaginary part of self-energy has a maximum near 
the Fermi energy (not shown); its position shifts toward 
Ep and its value approaches zero as the temperature 
is lowered. This suggests a Fermi-liquid like behavior 
for the model. Furthermore, the real part of the effec- 
tive self-energy for the d electrons shows a large nega- 
tive slope near the Fermi energy which denotes a large 
effective mass for the d electrons. We can use the self- 
energy to calculate the optical conductivity and the dc 
resistivity of the model (the details of calculations will 
appear in a future publication) . At temperatures of the 
order of Tq 1 the resistivity shows a log-linear behav- 
ior adjacent to a maximum, just like what is expected 
from a Kondo-type resistivity for a dilute system. As the 
temperature is lowered toward To the resistivity flattens. 
Below To the resistivity decreases quickly towards zero 
indicating the onset of coherence and the emergence of a 
Fermi-liquid (metallic) behavior. 

Interpretation Some of these results are consistent 
with a simple band-formation picture. When V = 
the available electronic states consist of a d band and 
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two (doubly degenerate) local / levels separated by U. 
When V is turned on, a new resonant state forms near the 
Fermi surface. Furthermore, the original d band mixes 
with the local / levels and the resonant level, giving rise 
to a renormalized band which has / character near the 
renormalized / level energies and has d character far from 
them as we see in Fig. 3. The Kondo states near Ep 
have mostly / character indicating that the / electrons 
themselves are involved in screening the local moments 
through hybridization with the d band. This is like the 
situation in a single-band Hubbard model, where the elec- 
trons within the band are responsible for screening the 
"local moments" in that band. 

We can understand the emergence of the two energy 
scales (Tq IAM and To) and the protracted screening of 
the moments in the PAM using the arguments of Nozieres 
0. He argued that since the screening cloud of a local 
magnetic moment involves conduction electrons within 
rpSiAM jrp p Q f -j-^e Fermi surface, only a fraction of the 
moments n e ff ~ N,i(0)Tq IAM may be screened by the 
conventional Kondo effect. He then proposed that the 
screened and unscreened sites may be mapped onto par- 
ticles and holes of an effective Hubbard model with local 
Coulomb repulsion U. The screening clouds hop from site 
to site and effectively screen all the moments in a dynam- 
ical fashion. The hopping constant of this effective model 
is strongly suppressed relative to t* by the overlap of the 
screened and unscreened states. Thus the Kondo scale 
of the effective model becomes much less than Tq 1 . 
This leads to the protracted screening behavior which is 
a crossover between the two regimes of Kondo screening 
at the higher scale T^ IAM and coherent screening at the 
lower scale Tq. In this argument the two scales depend 
on n e f f and hence there can be no universal relation be- 
tween them. 

Comparison with experiment For the purpose of com- 
parison with experiment we focus our attention on the 
data for CeBei3 compound because it has a large enough 
Kondo temperature that current high resolution appa- 
ratus should be able to probe its Kondo resonant state. 
CeBei3 is a mixed-valent heavy fermion compound with 
a small number of conduction electrons per site and a 
rather large Kondo temperature (Tr- ~ 400-ftT) so the 
PAM with the parameters we chose for n<j = 0.4 is 
suitable for comparison with experimental observations 
on this material (T w 320K if t* w 2eV). A series 
of photoemission experiments on single crystalline sam- 
ples of Ce compounds including CeBei3 have suggested 
a temperature dependence for the Kondo peak which is 
much weaker than what one expects from single- impurity 
model calculations ||, this is consistent with Fig. 1. 

Photoemission spectra obtained by choosing various 
incident angles on CeBei3 crystals have shown a disper- 
sive structure near the Fermi surface which is persistent 
up to temperatures of the order of 102V P5|- As we 
showed in Fig. 1 the PAM predicts a Kondo peak at these 



higher temperatures whereas the SI AM does not. Our 
results also show the dispersive nature of the quasiparti- 
cles as shown in Fig. 3. The photoemission observations 
show that the quasiparticle peak intensity near Ep is a 
minimum at the zone center and tends to develop as k de- 
viates from the zone center and reaches its highest value 
at the zone corner. This is consistent with what we see 
in Fig. 2. 

Conclusions We find that in the metallic regime 
where < 0.8 and rif sw 1: i) The PAM predicts a 
Kondo peak which has a weaker temperature dependence 
and persists up to much higher temperatures (in units 
of the Kondo scale) compared to the predictions of the 
SIAM. ii) Both the intensity and the spectral weight of 
the Kondo peak are larger in the PAM than in the SIAM. 
hi) The Kondo peak intensity does not scale like 1/V 2 
in the PAM. iv) The Kondo peak is dispersive in the 
PAM making a heavy quasiparticle band which crosses 
the Fermi surface. This can be interpreted as the heavy 
band emerging from an effective Hubbard model for the 
local screening clouds introduced by Nozieres. v) The 
Kondo peak below the Fermi energy starts to develop 
only as k deviates from the zone center, consistent with 
experiment, vi) There are two relevant energy scales 
for the PAM: the onset of screening scale Tq IAM , and 
the onset of coherence scale Tq which is strongly sup- 
pressed compared to the latter. This is also consistent 
with Nozieres arguments. 

These results provide a consistent interpretation of the 
Los Alamos photoemission experiments involving single- 
crystals, but not the results of the Bell Labs and Neucha- 
tel groups involving poly-crystals. Thus, a controversy 
persists. However, we note that near the insulating 
regime (rid ~ rif ~ 1) both the PAM Kondo scale (To) 
and screening rate are enhanced compared to those of 
the SIAM and Nozieres' arguments do not apply. In the 
narrow region where the system is doped away from the 
insulating state, both the screened local moment and the 
spectra show impurity-like temperature dependence. It 
is possible that this, or a more realistic model, including 
the effects of disorder, orbital degeneracy, or crystal field 
effects may provide a unifying interpretation of all of the 
photoemission spectra. 
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FIG. 1. Near-Fermi-energy (Ep = 0) structure of the 
/ density of states for the asymmetric PAM (left) and the 
symmetric SIAM (right). The model parameters are U = 1.5, 
V = 0.6, n d = 0.4 for the PAM and U = 2.75, V = 0.5, 
rid = 1.0 for the SIAM. The temperature dependence of the 
peak is universal for the SIAM and hence independent of the 
band fillings. Note that To refers to the two different Kondo 
scales of the PAM (left) and the SIAM (right). 
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FIG. 2. The PAM / and d spectral functions for different 
values of the momentum vector k for U = 1.5, V = 0.6, 
rif ~ 1.0, rid = 0.6 and T /To = 0.46. is the unrenormalized 
band energy. The Fermi energy is located at uj = 0. 
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FIG. 3. Band structure for the PAM. The model param- 
eters are the same as in Fig. 2. However the features persist 
up to T '/To ~ 10. The solid line shows the real part of the 
Green's functions poles vs. the unrenormalized band en- 
ergy. The symbols show the positions of the maxima in the / 
and d spectral functions. We characterize these peaks to be of 
/ character whenever Af(ek, Eh) > Ad{tk, Ek) or d character 
whenever Af(e h ,E k ) < A d (e k ,Ek) (cf. Fig. 2). 
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